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1.  Reasons 


To  increasing  the  safety  factor  of  the  space  vehicles,  Perspirable  Skin,  a  new 
designed  Thermal  Protection  System  (TPS)  was  proposed  by  bio-mimicking  the  process 
of  human  sweating.  Perspirable  Skin  is  arranged  in  a  ‘Peg  and  Hole’  fashion,  as  shown  in 
Fig.  1 .  The  RCC  panels  are  with  numerous  holes,  in  which  cores  (pegs)  made  of  another 
material  are  assembled.  The  core  material  has  smaller  Coefficient  of  Thermal  Expansion 
(CTE)  that  RCC.  As  temperature  increases,  a  gap  is  produced  between  the  RCC  skin 
and  cores  due  the  CTE  difference. 

Once  the  temperature  reaches  the 
working  temperature,  a  compressed 
gas  contained  in  tank  inside  the  space 
vehicle  blow  out  from  the  gap.  The 
cold  gas  flows  over  the  surface  and 
mixes  with  the  atmospheric  air  to 
compensate  the  frictional  heat 
between  the  vehicle  and  the 
atmospheric  air,  keeping  the  surface 
at  a  lower  temperature  than  that  of 
present  TPS.  After  the  vehicles 
landed  and  its  surface  temperature 
decreases,  the  gap  closes  again. 

Due  to  the  high  temperature  working  condition  and  low  CTE  requirement,  ceramic  is 
the  best  candidate  materials  to  manufacturing  the  core  of  Perspirable  Skin.  Powder 
Metallurgy  (P/M)  is  the  study  of  the  processing  of  powders.  It  is  the  most  common 
fabrication  method  for  most  ceramic  and  some  metal  engineering  components  [1], 
Compared  to  most  other  processing  methods,  P/M  is  cheaper.  Therefore,  P/M  was 
considered  to  be  the  best  candidate  method  to  fabricate  core  materials  for  Perspirable 
Skin. 

The  three  main  steps  of  P/M  to  convert  powders  into  useful  objects  are  1)  powder 
processing,  2)  forming  operations,  and  3)  sintering  [1],  The  powder  processing  may 
involve  milling  and  mixing  [2],  Many  methods  are  being  used  in  forming  operations  of 
powders  such  as  injection  molding,  slurry  techniques  and  compaction  [1],  Conventional 
uniaxial  powder  compaction,  which  is  performed  with  the  pressure  applied  along  one  axis 
using  hard  tooling  of  die  and  punches,  is  the  most  widely  used  method  to  compact 
powders.  Samples  after  compaction  are  called  “greens”,  which  are  the  precursor  for 
attaining  a  dense  body  by  sintering.  Sintering  describes  the  formation  of  bonds  between 
particles  close  to  their  melting  point.  It  can  occur  at  temperatures  below  the  melting  point 
by  solid-state  atomic  transport  phenomenon.  However,  some  instances  involve  the 
formation  of  a  liquid  phase  [1,  3],  On  a  microstructural  scale,  the  bonding  occurs  as 
cohesive  necks  grow  at  the  particle  contacts  [1-3], 


Fig.  1  Schematics  of  Perspirable  Skin  TPS 
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To  perform  the  compaction  and  sintering  as  well  as  monitor  the  sintering 
kinetics/densification  in  the  sintering,  a  Material  Test  System,  MTS  Insight  300  (MTS 
Systems  Corp.  USA)  and  a  Thermomechanical  Analysor  (TMA),  SETARAM  SETSYS 
Evoluation  95  (SETARAM  instrumentation,  France),  were  procured. 

2.  MTS  Insight  300 

2.1  Introduction 

MTS  Insight  systems  are  ideal  for 
performing  standard  tests  such  as  peel, 
tear,  sheer,  tensile,  compression,  and 
flex/bend  [4],  MTS  Insight  300  offers 
high-capacity  testing  capabilities  up  to  300 
kN.  As  shown  in  Fig. 2,  it  features  a 
floor-standing  frame  made  for  high-capacity 
strength  tests  such  as  tensile,  compression, 
and  characterization  of  higher-force 
materials.  The  motion  of  the  crosshead  is 
controlled  by  code  TestWorks™  Testing 
Software.  Offering  user-defined  crosshead 
speeds  as  well  as  advanced  control  modes 
such  as  load  and  strain,  MTS  Insight  300  can  analyze  material  properties  of  larger 
specimens.  Common  applications  include  metals,  ceramics,  high-strength  components, 
large  fasteners,  structural  materials,  and  composite  materials.  Operators  can  easily  run 
tests,  review  results,  and  report  and  archive  data  with  the  flexible  and  user  friendly 
TestWorks™  Testing  Software.  Table  1  shows  some  product  specifications  of  MTS 
Insight  300. 


Table  1 .  Product  Specifications  of  MTS  Insight  300  [4] 


Force  Capacity 

300  kN 

Crosshead  Travel 
Space 

1150  mm 

Weight 

1050  kg 

Height  x  Width  x 
Depth 

2440  x  1133x685 
mm3 

Frame  Stiffness 

300  kN 

Maximum  Test 
Speed 

500  mm/min 

Power 

6.0  kW 

Speed  Accuracy 

±  0.05%  of  full  speed 

Maximum  Test 

Speed 

500  mm/min 

Position  Accuracy 

0.01  mm 

Minimum  Test  Speed 

0.001  mm/min 

Position  Resolution 

0.001  mm 

Motor  and  Drive  System 

AC  Brushless,  3  ph  sinusoidal 

Fig. 2  MTS  Insight  300  system 
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2.2  Application  in  Perspirable  Skin  project 


MTS  Insight  300  is  used  to  compact 
various  ceramic  powders  to  forming  green 
samples.  Fig. 3  shows  the  positions  of  cross 
heads  and  the  uniaxiadie  set  in  the 
compaction  process.  The  detailed  compaction 
process  is  shown  in  Fig.  4.  The  loose,  free 
flowing  powder  is  introduced  into  the  die 
cavity  in  the  filling  step.  The  upper  surface  of 
the  powder  was  flattened.  Then  the  top  punch 
is  forced  into  the  die  to  compact  the  powder. 
After  compaction,  the  relative  density  of 
compact  is  around  50%  for  most  ceramic 
powders.  In  the  ejection  step,  the  forces 
acting  upon  the  punches  are  removed  and  the 
top  punch  is  used  to  eject  the  sample  from  the 
confines  of  the  die. 


Die  Set 


Cross 

heads 


Fig. 3  The  assembly  of  die  set  in  the 
MTS  Insight  300 


Except  fabricating  green  samples,  MTS  Insight  300  was  also  used  to  1)  investigate  the 
Powder  mixing  effect  on  the  compaction  capabilities  of  ceramic  powders;  2)  determine 
the  mechanical  properties  of  greens;  and  3)  investigate  the  density  distribution  inside  a 
green  compact. 


Filling  Flattening  Entry  of  Compaction  Ejection 

top  punch 


Fig. 4  The  process  of  die  set  compaction. 

2.3  Representative  experimental  results 

2.3.1  Powder  mixing  effect  on  the  compaction  capabilities  of  ceramic  powders  f51 

In  the  compaction  experiments,  by  recording  the  force,  displacement,  stress  and  strain, 
MTS  Insight  300  can  provide  us  the  relative  density-stress  curves  or  relative 
density-strain  curves.  Those  curves  were  compared  in  order  to  investigate  the  powder 
mixing  effect  on  the  compaction  capabilities  of  ceramic  powders  [5], 

In  the  series  of  experiments,  four  high-pure  and  undoped  alumina  (AI2O3)  powders 
with  submicron/micron  particle  sizes  were  investigated.  The  four  powders  were  TMDAR 
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(Taimei  Chemical  CO.,  LTD.,  Japan),  CR-15,  CR-6,  and  GE-1  (Baikowski  Ind.  Corp., 
U.S.A.).  All  possible  two-powder  combinations  with  the  four  powders  were  mixed  in 
various  mass  ratios,  and  in  total  six  powder  mixture  systems  were  produced.  The  mixture 
systems  include  TMDAR+CR-15,  TMDAR+CR-6,  TMDAR+GE-1,  CR-15+CR-6, 
CR-15+GE-1,  and  CR-6+GE-1.  In  addition  to  four  unmixed  powders,  each  system 
contains  nine  powder  mixtures,  in  which  the  proportions  of  one  powder  in  the  mixtures 
range  between  10%  and  90%  in  10%  intervals.  In  total,  fifty-four  powder  mixtures  were 
studied.  Six  grams  of  each  powder  mixture  poured  into  a  single-action  die  made  of  1144 
Stress  Proof  steel.  The  inner  diameter  of  the  die,  D/n,  is  22.19mm.  Before  pouring  the 
powders  into  the  die,  the  die  wall  was  lubricated  with  graphite  powder  (Panef  Corp., 
USA). 

The  top  surface  of  the  powder  was  flattened  and  then  the  top  punch  was  inserted  into 
the  die  until  it  touched  the  powder.  The  compaction  load  was  applied  on  the  top  punch  by 
MTS  Insight  300.  The  speed  was  lOmm/min  and  the  stress  resolution  is  0.01 3MPa. 
During  the  compaction,  only  the  top  punch  moved.  The  compressive  forces  applied  on 
the  top  (Ft)  and  bottom  punches  (Fb)  as  well  as  the  displacement  of  the  top  punch  (Zt) 
were  recorded  continuously  as  a  function  of  time  by  the  load  sensors.  At  first,  the 
pre-load  of  0.5MPa  was  applied  to  set  an  initial  state.  The  preload  is  necessary  to 
measure  the  compact  characteristics  accurately  before  starting  the  main  compacting 
process.  After  pre-loading,  the  initial  height  of  each  sample,  hi,  was  measured.  This  was 
done  by  measuring  the  total  height  of  the  two  punches  and  the  pre-loaded  compact 
inside  the  die,  then  subtracting  the  height  of  the  two  punches.  Considering  the  diameter 
of  each  compact  to  be  the  same  as  the  inner  diameter  of  the  die,  D/n,  the  initial  relative 
density,  R , ,  can  be  calculated  by  the  following  equation: 

p.  _  Pj  _  m/Vi  _  m 

Pi  pr  ^lxDPhl 

4 

where  m,  V,  and  p  represent  mass,  volume,  and  density,  respectively,  and  the 

subscripts  /  and  f  represent  initial  and  theoretically  fully  value.  The  theoretical 
fully-density  of  alumina  was  taken  as  3.975g/cm3. 

After  obtaining  the  initial  relative  density,  the  die  set  with  pre-loaded  compact  was  put 
back  into  the  load  frame,  and  then  the  load  of  80MPa  was  applied  to  compact  samples. 
The  final  green  samples  produced  a  diameter  of  22.19mm  and  their  heights  ranged  from 
7.05mm  to  9.1mm.  Several  green  samples  were  coated  with  nail  polish  and  then  their 
volumes  were  measured  in  water  by  Archimedes’  principle.  Those  values  were  then 
compared  to  the  corresponding  volumes  calculated  by  the  final  sample  heights,  hf,  and 
diameter,  D/n.  The  discrepancies  between  the  two  groups  of  values  were  minimial,  less 
than  1.1%.  Therefore,  only  the  calculated  volumes  were  used  to  calculate  the  densities. 

By  the  recorded  top  punch  displacement,  Zt,  and  the  initial  compact  height,  hi,  or  final 
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compact  height,  hf,  the  sample  height,  ht,  at  time  t  can  be  calculate, 


[  ht  -  hi  -  Zt  (2. a) 

\ht  =  hf  +  Zt  (2.b) 


The  ht  calculated  by  Eq.  (2. a)  and  (2.b)  for  any  compact  was  the  same.  Therefore,  two 
separate  measurements,  h ,  and  hf  for  each  compact  were  for  self-consistency. 

From  the  calculated  ht,  similar  to  calculating  the  initial  relative  density,  the  relative 
density  at  time  t,  defined  as Rt ,  can  be  found  as, 

n  _  pt  _  m/Vt  _  m  (3) 

4 

In  Eq.  (3),  the  sample  diameter  was  assumed  to  be  matching  the  inner  diameter  of  the 
die,  Din. 

The  corresponding  stress,  crt,  at  time  t  can  be  obtained  by  the  punch  diameter,  D/n, 
and  the  recorded  load,  Ft. 


_Ft  _ 

Gt  St  1 


Ft 


71  Din 


(4) 


where  St  is  the  cross-sectional  area  of  the  top  punch. 


The  relative  density-stress  curves  for  each  powder  mixture  group  were  then  drawn. 
Fig.  5. a  and  Fig.  5.b  show  the  representative  curves  for  the  TMDAR+CR-15  and 
CR-6+GE-1  powder  mixture  groups,  respectively.  For  clarity,  only  six  curves  were 
included  in  each  figure. 


Fig. 5  The  relative  density-stress  curves  by  the  compaction  for  a)  CR-15+TMDAR  and  b) 
CR-6+GE-1  [5], 

Under  a  constant  compaction  load  or  stress,  the  difference  of  the  relative  densities 
represents  the  change  in  the  compaction  capability  in  different  mixtures.  A  higher  relative 
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density  represents  improved  compaction  capability.  Even  the  absolute  relative  density 
values  continuously  increase  with  the  increase  in  the  compaction  load,  the  relationship 
between  the  relative  density  and  powder  proportions  remains  similar  for  each  powder 
mixture  group  as  observed  in  Fig.  5. 

Therefore,  one  can  conclude  that  the 
starting  powder  combination  determines 
the  compaction  capability  of  each  powder 
mixture.  In  addition,  the  initial  relative 
densities  dictated  by  the  starting  powder 
combination  can  be  used  to  quantify  and 
compare  the  compaction  capability  of 
each  powder  mixture. 

The  relationships  between  the  initial 
relative  density  and  the  powder 
proportion  are  shown  in  Fig. 6  for  all  six 
powder  mixture  groups.  For  each  curve, 
the  left  end  is  the  100%  of  the  fine 
powder  and  the  right  end  is  the  100%  of 
the  coarse  powder.  The  “fine”  or  “coarse” 
powder  is  defined  by  the  average  particle 
size  D  of  the  powder. 

Dge  -  i  >  Dcr  -  6  >  Dcr  - 15  >  Dtmdar 

Fig.  6  shows  two  kinds  of  relationships  between  the  initial  relative  density  and  the 
proportion  of  the  coarse  powder. 

1.  Linear  relationship:  TMDAR+CR-15  and  CR-15+CR-6  powder  mixture  groups 
show  this  characteristic.  The  maximum  relative  density  corresponds  to  100%  of 
one  of  the  powders. 

2.  Nearly  parabolic  relationship:  The  maximum  relative  density  of  some  powder 
groups  with  this  relationship  occurs  with  certain  powder  mixtures,  such  as 
TMDAR+CR-6,  CR-15+GE-1,  and  CR-6+GE-1.  This  is  similar  to  the  conclusion 
of  Westman  [6]  and  McGeary  [7],  However,  the  proportions  of  coarse  powder  in 
the  powder  mixture  that  attained  the  maximum  relative  density  are  between  40% 
and  60%.  These  values  differ  from  the  reported  30%  by  Westman  [6]  and 
McGeary  [7]  for  sand  and  metal  balls.  Another  powder  systems,  TMDAR+GE-1, 
has  the  maximum  relative  density  at  100%  TMDAR. 

In  the  study  of  ceramic  powders’  compaction,  Aydin  et  al.  [8-10]  defined  axial  true 
strain  s  t  as  the  following, 


Fig. 6  The  initial  relative  density  -  coarse 
powder  proportion  for  all  six  powder 
mixture  groups  after  0.5MPa  pre-load 


compaction  [5], 

(5) 


st- 


<) 


(6) 
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Fig.  7(a)  and  Fig.  7(b)  show  the  representative  relative  density-strain  curves  for  the 
TMDAR+CR-15  and  CR-6+GE-1  powder  mixture  groups,  respectively.  Comparing  Fig.  7 
and  Fig.  5,  the  powder  mixture  with  better  compaction  capability/density  in  Fig.  5  is  also 
with  a  higher  relative  density  under  any  strain  values  in  Fig.  7. 


Fig. 7  The  relative  density-strain  curves  by  the  compaction  for  a)  CR-15+TMDAR  and  b) 
CR-6+GE-1  [5], 

The  Mircocal  Origin  and  Wolfram  Mathematica  software  were  employed  to  calculate 
the  trend  functions  of  the  relative  density-strain  curves.  The  relationship  exactly  fits  into 
the  following  equation  for  each  powder  mixture, 

Rt  =  KQst  (7) 

where  k  has  the  same  unit  as  density,  g/cm3.  Table  2  displays  several  k  value  for  the 
TMDAR+CR-15  powder  mixture  group. 

Table  2  The  coefficient  k  in  Equation  (3.7)  for  representative  TMDAR+CR-15  powder 
mixture 


CR-15  0%  10%  50%  70%  90%  100% 

(g/cm3)  1.4228  1.3843  1.1887  1.0688  0.9622  0.908 


Eq.  (7)  shows  that  a  powder  mixture  with  a  higher  density/compaction  capability  has  a 
larger  k  value.  Compared  the  k  values  in  Table  2  to  the  corresponding  initial  relative 
density  value  in  Fig. 6,  they  are  the  same.  This  indicates  the  physical  meaning  of  k  in  Eq. 
(7)  is  the  “initial  relative  density”.  From  the  relative  density  and  strain  definitions  in  Eq.  (3) 
and  (6)  as  well  as  Eq.  (1),  the  conclusion  can  be  verified. 

Therefore,  the  method  to  compare  compaction  capability  of  various  powder  mixtures 
from  their  relative  density-strain  curves  only  needs  to  consider  the  initial  relative  density 
value.  In  other  words,  from  either  the  relative  density-stress  or  the  relative  density-strain 


curves,  one  can  use  Fig.  6  to  find  the  difference  in  the  compaction  capabilities  of  powder 
mixtures. 


By  further  analysis,  we  found  the  shape 
of  the  initial  relative  density  versus 

proportion  of  coarse  powder  curve 

depends  on  the  ratio  of  average  particle 
size,  and  the  difference  in  the  initial  density 
of  two  powders.  If  the  average  particle  size 
ratio  is  less  than  3.6,  or  the  average 
particle  size  ratio  is  larger  than  3.6  but  the 
initial  relative  density  different  is  higher 
than  10%,  the  curve  is  linear.  If  the 

average  particle  size  ratio  is  larger  than  3.6 
and  the  initial  relative  density  difference  is 
smaller  than  10%,  the  curve  is  parabolic. 
The  triangle  enclosed  by  the  line 

intersecting  the  relative  density  of  the  fine 
powder  and  the  maximum  relative  density 
of  the  coarse  powder,  the  line  intersecting 


0  20  40  60  80  100 

Proportion  of  coarse  powder  (%) 


Fig. 8  Triangle  boundaries  for  the  initial 
relative  density-coarse  powder  proportion 
curves  of  TMDAR+CR-6,  CR-6+GE-1 
powder  groups  [5], 


the  relative  density  of  the  coarse  powder  and  the  maximum  relative  density  of  the  fine 


powder  as  well  as  the  line  connecting  the  relative  density  of  the  fine  and  coarse  powder 


form  the  boundary  of  the  initial  relative  density-proportion  of  coarse  powder  curve,  as 


shown  in  Fig. 8  [5], 


2.3.2  Determination  of  the  mechanical  properties  of  greens 


The  stress-strain  relationship  in  both 
loading  and  unloading  parts  can  be 
recorded  by  MTS  Insight  300,  as  shown  in 
Fig. 9.  The  elastic  unloading  part  of  the 
stress-strain  curves  for  the  powders  is  linear. 
The  slope  of  the  elastic  unloading  part  ( m ), 
the  bulk  modulus  (K),  the  shear  modulus 
(G),  the  Young’s  modulus  (E)  and  the 
Poisson’s  ratio  (o)  have  relationship  as 
listed  in  Eqs.  8-10  [11], 

m=K  +  -G  (8) 
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Strain 

Fig. 9  Stress-strain  curve  of  TMDAR  in 


both  loading  and  unloading  parts 


9  GK 
3  K  +  G 


(9) 


9 


3K  -2G 
V  ~  2(3 K  +  G) 


(10) 


The  TMA  unit  offers  a  facility  to  measure  Young’s  modulus  of  materials  (E)  using  the 
three-point  bending  method,  which  will  be  discussed  in  Section////.  The  Young’s  modulus 
of  CR-15  compact  by  80MPa  is  1031 .3  MPa,  and  that  of  TMDAR  is  2155  MPa.  Therefore, 
by  measured  Young’s  modulus  and  the  calculated  slope  value,  the  other  moduli  and 
Poisson’s  ratio  can  be  obtained,  as  listed  in  Table  3. 


Table  3  Material  properties  for  several  alumina  compacts 


Powder 

TMDAR 

CR-15 

Young’s  modulus  (MPa) 

2155.1 

1031.3 

Unloading  slope  (MPa) 

2498.3 

1250.4 

Shear  modulus  (MPa) 

876.1 

410.8 

Bulk  modulus  (MPa) 

1329.9 

702.3 

Poisson’s  ratio 

0.23 

0.255 

2.3.3  Density  distribution  inside  green  compacts 

The  densities  of  a  green  compact  are  not  even  everywhere,  but  with  a  density 
distribution  [8,  9],  We  used  colored  powder  as  tracer  to  investigate  the  density  distribution 
insider  green  compacts. 

In  Section  2.3.1,  we  already  discovered  that  powder  mixing  changes  the  compaction 
capability.  Therefore,  we  dyed  several  powder  mixtures  by  stamp-ink  (2000  plus,  Red, 
COSCO,  Germany).  The  mixtures  then  were  put  into  a  furnace  (Carbolite-HTF1700,  UK) 
for  4  hours  at  70°C  to  dry  out.  The  remained  solids  were  put  into  bottles  with  12mm 
diameter  alumina  mixing  media.  The  bottles  were  on  a  jar  mill  (U.S.  Stoneware  764AVM, 
U.S.A.)  for  48  hours  to  break  the  solids.  Then  the  powders  were  sifted  to  get  colored 
powder  with  size  less  than  0.3mm. 

The  colored  powders  were  compacted  to  0.5MPa  to  get  the  initial  densities  as 
mentioned  in  Section  2.3.1.  Then  the  values  were  compared  with  five  different  main 
powder  (mixture):  TMDAR,  50/50  TMDAR/CR15,  CR15,  50/50  CR15/GE1,  and  GE1, 
respectively.  The  colored  powder  with  same  compaction  capability  to  a  powder  (mixture) 
will  be  used  as  its  compaction  indicator. 

To  investigate  the  density  distribution  insider  a  green  compact,  the  height  of  the  green 
sample  cannot  be  two  small.  Otherwise,  the  density  distribution  will  be  in  a  small  range, 
hard  to  identify.  So  in  this  study,  we  keep  the  height/diameter  ratio  of  each  green  sample 
around  1.1-1. 2.  The  diameter  of  the  samples  is  19.04  mm.  In  the  experiments,  we  kept 
the  thickness  of  each  main  powder  layer  are  identical,  also  each  colored  powder  layer. 
And  the  thickness  of  colored  powder  layer  is  around  1/15  of  the  main  powder  layer. 
Based  on  these  requirements,  the  powder  mass  of  each  layer  can  be  calculated.  The 
main  powder  is  1.5g  for  each  layer  and  the  colored  powder  mass  is  between  0.1 2g  for 
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each  layer. 

Similar  to  the  compaction  procedure  in  Section  2.3.1,  the  first  layer  powder,  which 
always  is  main  powder,  was  pour  into  the  die,  which  was  lubricated  with  graphite  powder 
(Panef  Corp.,  USA).  Then  a  pre-load  of  0.5MPa  was  applied  to  compact  the  first  layer  to 
obtain  enough  stiffness.  After  the  pre-compaction,  the  die  set  was  moved  from  the  MTS 
Insight  300.  The  layer  thickness,  measured  by  subtracting  he  height  of  the  two  punches 
form  the  total  height  of  the  two  punches  and  the  pre-loaded  compact  inside  the  die,  was 
recorded.  The  top  punch  was  then  dragged  out  of  the  die  slowly  and  then  the  second 
layer  powder  (colored  powder)  was  pour  into  the  die. 

Flatting  the  surface  of  the  new  layer  powder  and 
putting  the  top  punch  back  to  the  die,  the  0.5  MPa 
pre-load  was  applied  on  the  both  layer  again.  Repeat 
the  procedure  till  the  last  layer  was  pre-loaded. 

Finally  the  load  of  80MPa  was  applied  to  compact 
the  whole  multi-layered  samples. 

The  green  samples  were  then  took  out  the  die 
and  polished  carefully  on  the  lateral  surface  till  to  the 
middle  section.  Then  each  sample  was  placed  on  a 
sample  holder  for  photo  taking.  The  photos  were 
took  by  a  dissection  microscope  (Wild  M5A,  Wild 
Heerbrugg  Ltd,  Switzerland).  The  representative 
picture  of  the  middle  section  is  shown  in  Fig.  10.  Due 
to  the  axial  symmetric  characteristics  of  the  density 
distribution  in  the  compact,  the  broken  parts  in  the 
left  edge  shown  in  Fig.  10  do  not  influence  the 
calculation  of  the  density  distribution  contour. 

Based  on  Eqs  6  and  7,  the  relative  density  pt  can  be  calculated  as: 


Fig.  10  The  photo  of  the  middle 
section  of  a  TMDAR  sample.  The 
bottom  side  in  the  picture  is  also 
the  side  contact  to  the  bottom 
punch  in  the  compaction. 


pt  =  ptQ 


ln# 


(11) 


where  p,  is  the  initial  density  after  0.5MPa,  which  already  is  presented  in  Fig. 6  for  each 
main  powder.  The  final  height  of  a  position,  ht,  is  obtained  by  analyzing  the  photo  by 
software  Image  G.  Since  each  indicator  layer  is  curved  after  compaction  as  shown  in 
Fig.  10,  the  ht  varies  by  positions.  Therefore,  different  position  is  with  different  density. 

The  density  distribution  contour  was  plotted  by  MatLab,  Fig.  11  shows  the  contour  of 
TMDAR,  in  which  the  0  position  is  corresponding  to  the  center  axial  in  the  Fig.  10.  All 
other  four  powder  compacts  feature  similar  density  distribution  contour  as  shown  in 
Fig.11.  But  the  density  distribution  range  varies  as  different  powders  as  listed  in  Table  4. 
Similar  to  the  average  density  relationship  in  Fig. 6,  powders  with  better  compaction 
capability  have  both  larger  lowest  density  and  highest  density.  From  Fig.  11,  one  can  see 
that  the  high-density  regions  exist  at  the  near  the  bottom  central  region  and  on  the  central 
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axis.  And  the  low-density  regions  are  at  the  bottom  edges  and  near  the  top  center.  Our 
results  are  similar  to  Aydin  et  al  [8,  9]  and  others’  report  [13-15], 


Fig.  11  The  density  distribution  of  TMDAR  green  sample,  zero  position  is  the  radial 
symmetric  axel  in  Fig.  10. 


Table  4  The  lowest  and  highest  density  in  various  samples 


Powder 

Lowest  density  (g/cm3) 

Highest  density  (g/cm3) 

TMDAR 

1.618 

2.462 

CR-15 

1.199 

1.662 

GE-1 

1.203 

1.671 

50/50  TMDAR/CR15 

1.409 

2.064 

50/50  CR15/GE1 

1.414 

1.867 

3.  TMA  -  SETARAM  SETSYS  Evoluation  95 

3.1  Introduction 

SETARAM  SETSYS  Evoluation  95  is  designed  by  SETARAM  instrumentation,  France 
(http://www.setaram.com/).  This  machine  measures  the  deformation  of  a  sample  under 
non-oscillating  stress  against  time  or  temperature  with  prescribed  ramp/soak  path.  The 
testing  temperature  range  of  SETSYS  Evolution  95  is  from  room  temperature  to  1750°C 
[16].  Fig.  12  shows  the  front  view  of  this  analyzer. 

When  the  head  cover  open,  the  sample  can  be  put  in  sample  chamber,  which  is  an 
alumina  tube.  A  probe  connecting  to  a  displacement  transducer  is  arranged  in  the  sample 
chamber.  The  position  of  the  probe  can  be  adjusted  so  its  tip  can  touch  the  top  surface  of 
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the  sample,  as  shown  in  Fig.  13.  For  the 
sintering  kinetics  and  CTE  test,  a  5-gram  load  is 
applied  to  the  probe  to  ensure  the  contact  with  a 
sample  throughout  the  entire  experiment.  It  has 
been  verified  that  this  small  force  does  not 
influence  the  testing  results  [17].  For  the 
Young’s  modulus  measurement,  the  force  can 
be  as  large  as  200  grams.  The  displacement 
transducer  equipped  in  the  SETSYS  Evolution 
95  is  characterized  by  its  robustness  and  high 
accuracy;  it  can  detect  dimension  changes  as 
small  as  0.01  micron.  The  transducer  uses  an 
electromagnetic  system  for  automatic  control  of 
the  force  applied  to  the  sample,  between  0.01 
and  1.5  N.  The  force  can  be  increased  by 
adding  weights  (up  to  200  grams)  on  a  top  plate. 
Transducer  calibration  and  force  control  are 
managed  automatically  by  computer  [16]. 

The  SETSYS  Evoluation  95  has  a  cylindrical 
furnace.  The  heating  element  is  made  up  of 
graphite  rube  and  fitted  in  the  centerline  of  the 
furnace.  The  thermocouple  is  composed  of  Pt/Pt 
high  as  1750°C  [16]. 


Fig.  12.  The  front  view  of  the  SETARAM 
SETSYS  Evoluation  95  (1-head  cover, 
2-furnace  side  watch,  3-mode  changing 
switch,  4-open/close  head  switch, 
5-side  cover,  6-proctective  gas  inlet  and 
outlet  tube). 

-Rh  and  can  withstand  temperatures  as 


3.2  Application  in  Perspirable  Skin  project 

Except  as  a  regular  furnace  to  sinter  samples  with 
oxidization-protection  environment,  with  its  software, 
the  SETSYS  Evolution  95  can  be  used  as  a 
dilatometer  to  test  many  material  properties,  such  as 
the  thermal  expansion  coefficients  (CTE),  the  softening 
temperature,  the  glass  transition  point  and  so  on.  The 
SETSYS  Evolution  95  is  also  particularly  suitable  for 
measuring  controlled-rate  sintering.  By  measuring  the 
displacement  of  one  dimension  and  performing 
calculations,  this  TMA  can  give  the  shrinkage 
rate-sintering  temperature  curve  or  densification 
rate-relative  density  curve.  By  a  specially  designed 
probe  and  support,  the  SETSYS  Evolution  95  also  can 
perform  3-point  flexure  measurements,  which  would 
give  the  temperature  dependence  of  the  Young’s 
modulus  of  the  testing  material.  Fig.13  Sample  position  inside 

the  TMA 
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3.3  Representative  experimental  results 

3.3.1  Sintering  kinetics  curves  [18,  191 


Two  types  of  curves  are  commonly  used  to  describe  the  sintering  kinetics:  the 
shrinkage  rate-sintering  temperature  curve  and  the  densification  rate-relative  density 
curve.  For  the  latter  one,  the  relative  density  D  is  the  ratio  of  the  mass  density  of  the 
sample  to  the  theoretical  density  of  the  corresponding  powder  mixture.  And  the 
densification  rate  is  defined  as: 


s\T) 


_  dp  _  D(T ) 
pdt  D(T) 


(12) 


We  investigated  the  sintering  kinetics  for  a  series  of  alumina  and  zirconia  powder 
mixtures  by  TMDAR,  CR-15  (alumina),  as  well  as  TZ3YS  and  CERAC-2003  (zirconia). 
The  proportion  of  TMDAR  in  alumina  system  and  CERAC-2003  in  zirconia  system 
ranges  from  0%  to  100%  in  10%  intervals,  resulting  in  a  total  of  22  powder  mixture 
samples  in  this  study. 

From  the  TMA  testing  results,  we  calculated  and  plotted  densification  rate-relative 
density  curves  to  describe  the  sintering  kinetics  of  each  powder  mixture.  Since  the  TMA 
can  measure  the  length  change  in  only  one  direction,  to  calculate  the  densification  rate 
based  on  the  results  of  TMA,  Lance  et  al.  [20]  defined  an  anisotropic  shrinkage  factor  cr: 


{(pf-(p0)L0 

a  =  — ro>  (13) 

(^7  ~  L0  )(po 

where  Lo  and  Lf  are  the  initial  and  final  height  of  the  specimen,  and  cpo  and  cpf  are  the 
initial  and  final  mean  diameters  of  a  specimen.  For  each  of  our  experimental  results,  a 
was  almost  1  (>0.991 ).  It  is  therefore  reasonable  to  consider  the  sintering  process  for  our 
samples  to  be  a  self-similar  (isotropic)  one.  Thus,  with  the  1:1  height  to  diameter  ratio, 
the  temperature  dependent  relative  density  can  be  calculated  from  the  following 
expression: 


D(T)  = 


i+7_7 

L„ 


1  + 


A L(T) 


Df 


L, 


L0  +  A L(T) 


Dt 


(14) 


where  Df  is  the  relative  density  of  the  fully-sintered  sample,  AL(T)  is  the  displacement  at 
temperature  T,  so  the  L0+AL(T)  represents  the  height  of  the  sample  at  temperature  T. 
Then  it  is  relatively  easy  to  derive  the  expression  of  densification  rate: 


,sm_D(T)_  -A L{T) 

K  ’  D(T )  20[Lo  +  AL(T)] 


(15) 
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where  the  coefficient  1/20  (3/60)  is  a  result  of  the  conversion  to  SI  units,  since  the  TMA 
counts  time  in  minutes  instead  of  seconds. 

Fig.  14  shows  the  representative 
densification  rate-relative  density  curves 
observed  in  the  alumina  powder  mixtures. 

Because  of  the  similarity  of  the  curves  for  1 ) 

10%-30%  TMDAR  powder  mixtures,  2) 

60%-90%  TMDAR  powder  mixtures,  only  six 
densification  curves  including  0%  and  100% 

TMDAR  powders  are  shown  to  enhance  the 
clarity  of  the  plots. 

As  shown  in  Fig. 14,  increasing  the 
proportion  of  TMDAR  enhances  both  the 
green  density  and  the  final  density  of  the 
samples.  With  the  exception  of  the  pure 
CR-15  sample,  all  of  the  specimens  are 
nearly  fully-sintered  (D  >  95%).  Due  to  the 
high  sintering  capability  of  TMDAR,  the  samples  containing  a  higher  proportion  TMDAR 
have  higher  fully-sintered  densities. 

In  the  densification  rate-relative  density  curve,  one  point  of  critical  importance  is  the 
maximum  densification  rate  point.  Prior  to  reaching  this  point,  there  is  no  grain  growth 
and  the  pores  among  the  grains  are  parts  of  interconnected  networks.  These  networks 
begin  to  collapse  into  isolated  pores  after  the  powder  compact  reaches  its  maximum 
densification  rate,  and  rapid  grain  growth  begins  to  occur  [21].  From  Fig.  4.8,  the 
maximum  densification  rate  for  the  alumina  powder  mixtures  occurs  within  the  relative 
density  range  of  70-84%,  consistent  with  the  results  of  Lance  et  al.  [20]  and  Lange  [22], 
Another  point  of  interest  is  that  the  increase  in  the  proportion  of  TMDAR  also  increases 
the  relative  density  at  which  the  maximum  densification  rate  occurs.  This  change  is  also  a 
result  of  the  increased  green  density  of  the  samples  caused  by  increased  TMDAR 
proportion. 

The  ascending  and  descending  parts  of  a  densification  curve  are  the  segments  before 
and  after  the  maximum  densification  point,  respectively.  Fig.  14  shows  that  the  powder 
mixtures  with  higher  CR-15  proportions  have  higher  densification  rates  in  the  ascending 
section  and  lower  rates  in  the  descending  section  of  the  curve.  Ting  et  al.  [21]  concluded 
by  theoretical  calculation  that  “Prior  to  the  occurrence  of  grain  growth,  a  maximum 
densification  rate  exists  as  the  starting  particle  size  distribution  width  increases,  . . .  During 
grain  growth,  the  model  predicts  decreasing  densification  rate  with  increasing  starting 
particle  size  distribution  width."  Our  experimental  observations  match  this  conclusion 
since  the  increase  in  the  proportion  of  CR-15  increases  the  particle  size  distribution 
width. 

Depending  on  the  powder  mixture  system,  the  ascending  section  for  the  densification 


Relative  density 


Fig.  14  Densification  rate  as  a  function  of 
relative  density  for  TMDAR+CR-15 
alumina  powder  mixtures.  [18] 
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rate-relative  density  curves  can  be  closely  approximated  by  either  one  or  two  straight 
lines.  The  100%  TMDAR,  agglomerate-free  sample,  exhibits  a  single  linear  segment  with 
one  slope.  In  contrast,  systems  containing  CR-15  demonstrate  two  linear  segments  with 
distinct  slopes.  Lance  et  al.  [20]  obtained  similar  results  by  stating  “Agglomerates  free 
alumina  powder  shows  a  single  slope  linear  segment  and  powders  with  agglomerates 
show  two  linear  segments 

Each  powder  mixture  containing  CR-15  features  a  plateau  in  the  ascending  portion  of 
the  curve  in  Fig.  14.  Even  though  the  relative  density  and  densification  rate  value  are 
different  for  each  powder  mixture,  the  sintering  temperature  at  which  the  plateau  occurs 
is  within  the  same  temperature  range  (1150-1200°C).  This  plateau  is  a  result  of  y-  to 
a-phase  transition  of  alumina  [23],  since  10%  y-phase  alumina  is  present  in  CR-15  (Table 
2.1).  TMDAR  powder  contains  only  a-phase  alumina  and  no  plateau  in  the  pure  TMDAR 
densification  rate-relative  density  curve. 

A  plateau  in  descending  section  of  a  densification  rate-relative  density  curve  signifies 
the  occurrence  of  abnormal  grain  growth  [20],  It  can  be  seen  that  no  plateaus  exist  in  that 
region  for  any  curves  in  Fig.  4.8.  This  holds  true  even  for  a  powder  system  with 
agglomerates  (such  as  CR-15),  which  is  known  to  cause  abnormal  grain  growth. 


Fig.  15  shows  the  seven  representative 
curves  obtained  for  our  zirconia  powder 
system.  The  10%  and  20%  CERAC-2003 
systems  as  well  as  all  the  powder  mixture 
systems  between  50%  and  80% 

CERAC-2003  exhibit  very  similar 
densification  rate-relative  density  curves. 

The  densification  rate-relative  density  curve 
for  each  group  is  presented  in  Fig.  15  as  a 
single  curve  to  allow  for  clarity. 

Fig.  15  shows  that  the  green  density  of 
zirconia  samples  increases  when  the 
proportion  of  CERAC-2003  is  between  0% 
and  40%,  and  begins  to  decrease  beyond 
40%.  This  is  because  the  small  TZ3YS 
particles  fill  in  the  interstitial  spaces  among  the  larger  CERAC-2003  particles  as  well  as 
the  powders  of  the  zirconia  system  feature  very  similar  compactibility.  The  similarity  in  the 
sinteribilities  of  CERAC-2003  and  TZ3YS  results  in  a  similar  final  density  for  various 
zirconia  powder  mixtures. 


0.6  0.7  0.8 

Relative  density 

Fig.  15  Densification  rate  as  a  function  of 
relative  density  for  TZ3YS+CERAC-2003 
zirconia  system.  [18] 


The  maximum  densification  rate  for  each  of  our  zirconia  powder  mixtures  occurs 
within  the  relative  density  range  of  74-82%.  This  result  is  consistent  with  the  report  of 
Trunec  et  al.  [24],  The  higher  the  green  density  of  a  powder  mixture,  the  higher  the 
relative  density  at  which  the  maximum  densification  rate  occurs.  This  is  the  same  result 
found  with  our  alumina  powder  mixtures.  The  densification  rate  in  ascent  and  descent 
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part  changes  as  the  proportion  of  CERAC-2003,  which  also  fit  to  the  Ting’s  calculation 

[21]. 

3.3.2  Phenomenological  constitutive  models  for  sintering  curves  [18,  19] 


In  the  last  decades,  a  lot  of  research  has  been  devoted  to  predict  the  densification 
kinetics  of  powder  compaction  [25-27],  Some  of  this  research  used  the 
phenomenological  model.  The  phenomenological  model  consists  of  fitting  analytical 
expressions  directly  from  the  results  of  sintering  experiments.  Thus,  it  provides  a 
constitutive  equation  describing  as  precisely  as  possible  the  real  behavior  of  the  material 
without  regarding  the  physical  meaning  of  its  parameters. 


An  example  of  the  phenomenological 
models  is  presented  by  Hsueh  [27]: 

ss=Q.(T)[Dn(T)-Df  (16) 

where  D~(T)  is  the  theoretically  possible 
full-sintered  density  at  temperature  T. 

To  fit  the  densification  rate  -  relative 
density  curves  to  Eq.  (16),  part-isothermal 
heating  cycle  sintering  of  samples  is 
needed.  The  part-isothermal  heating  cycle 
sintering  was  first  utilized  by  Dorn  in  1957 
[28],  The  part-isothermal  heating  cycle 
used  in  ours  study  is  shown  in  Fig.  16. 

Fig.  17  shows  the  part-isothermal 
sintering  test  result  for  100%  TMDAR 
powder.  Each  “descending  part”  is  from  the 
isothermal  sintering.  If  one  extends  each 
descending  part  and  intersects  the  extension 
with  the  relative  density  axis,  the  intersection 
point  gives  the  D«,  for  the  temperature  at 
which  the  isothermal  sintering  performed. 

After  obtaining  a  series  of  values  of  Dco(T), 
the  following  steps  are  used  to  calculate  the 
phenomenological  constitutive  model  for 
each  powder  mixture  in  Figs.  4.8  and  4.9: 

1 .  Taking  the  natural  logarithm  for  the 
both  sides  of  Eq.  (16),  a  linear 
relationship  can  be  found: 
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Fig.  16  The  part-isothermal  ramp/soak 
sintering  cycle  used  in  this  study  to 
calculate  the  phenomenological  constitutive 
models.  [18] 
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Fig.  17  The  densification  rate  -  relative 
density  relationship  of  the  100%  TMDAR 
powder  under  the  part-isothermal 
sintering  path  described  in  Fig.  16.  [18] 
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\n(ss )  =  n  ln([Z)oo(r)  -£>])  +  ln(Q  (7)) 


(17) 


2.  Calculating  D«,(T)  -  D  for  the  relative  density  data  of  each  isothermal  sintering 
testing  part  and  then  taking  the  natural  logarithm  for  the  densification  rate  data 
and  the  D°°(7)  -  D. 

3.  Using  the  linear  function  to  fit  the  terms  ln(<£“v)and  ln(D°°(7")  -  D),  a  series  of 

linear  equations,  each  of  which  relates  to  a  specific  isothermal  sintering 
temperature  are  obtained.  Table  5  shows  the  linear  equations  calculated  for  pure 
TMDAR  powder. 


Table  5.  The  linear  equations  calculated  for  pure  TMDAR  powder  [19] 


Isothermal  temperature  (°C) 

Linear  fitting  equations 

1020 

\n(ss )  =  1 .63 8  lnr0.6 1 1 8  -  D\ - 1 .6439 

1080 

ln(£-s)  =  1.6521n[0.6802-D]  +  0.3165 

1140 

ln(^)  =  1 .648  ln[0.792-  D]-  0.0657 

1200 

ln(^)  =  1 .655  ln[0.973  -  D] -  0.5 1 73 

1260 

ln(^)  =  1 ,6351n[0.995  -D]- 1 .059 

4.  The  coefficients  in  front  of  ln(D°°(7")  -  D)  give  the  value  of  n  in  Eq.(1 7).  From  Table 
5  we  can  see  that  the  n  value  is  nearly  constant  for  a  powder  system  under 
different  temperature.  The  n  for  TMDAR  is  around  1 .65. 

5.  The  entries  in  Table  5  give  a  series  value  of  ln(Q(7")).  Thus  we  can  get  a 
relationship  between  ln(Q(T))  and  temperature  T,  then  use  a  polynomial  to  relate 
Q(7)  and  T.  For  pure  TMDAR,  the  Q(7")  can  be  expressed  as: 

Q(r)  =  -1.3799xl0-nr4+6.0705xl(T6 7 8r3 

(18) 

-8.719  x  l(T5r2  +  0.037079 T  +  5.2474 

6.  A  series  D«,{T)  values  at  different  temperature  T are  also  obtained,  we  can  fit  D«>(T) 
into  a  hyperbolic  tangent  function.  For  the  pure  TMDAR  powder,  the  expression  is: 

Doo(T)  =  (0.2224  tanh[8.3785  x  10“37  -  9.4498]}  +  0.7776  (19) 

7.  Combining  the  result  of  steps  4-6,  The  final  phenomenological  constitutive  model 

for  each  powder  mixture  is  obtained. 
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The  calculated  phenomenological  constitutive  function  for  100%  TMDAR: 

=  [-1.138x1 0-n(T)4  +  0.503408  x  1 0"7(r)3 
-  0. 72765  xl0~4(7")2  +0.03 19064(7") +  3. 5792]  (20) 

x{  (0.2224  tanh[8. 3785  xl0~3(7"  +  40)  -  9.4498]  }  +  0. 7776 -29}1 5 

Fig.  18  shows  the  comparison  of  the  experimental  measurement  and  the  calculated 
result  by  Eq.  (20)  for  both  the  non-isothermal  and  part-isothermal  heating  cycles.  The 
correlation  between  the  calculated  and  experimental  results  remains  very  high  for  both 
two  figures.  The  errors  between  the  experimental  data  and  the  calculated  results  are  less 
than  ±6%. 


Fig.  18  Comparison  of  the  experimental  and  calculated  results  for  100%  TMDAR  under  a) 
non-isothermal  and  b)  part-isothermal  heating  circle  [18] 


The  Hsueh’s  formula  also  can  be 
applied  to  zirconia  powder  systems.  The 
phenomenological  constitutive  equations 
for  all  zirconia  powder  mixtures  are 
obtained.  Fig.  19  shows  the  comparison  of 
the  experimental  and  the  calculated  results 
by  Eq.  (4.10)  during  the  non-isothermal 
heating  cycle.  The  errors  between  the 
experimental  data  and  the  calculated 
results  are  less  than  ±4%. 

The  equation  to  describe  the  sintering 
kinetics  of  pure  TZ3YS  is: 


Fig.  19  100%  TZ3YS  sintering  [18] 


ss  =  [1 .00465  x  1 0 -\T  - 1 50)4  +  5 .83063  x  1 0 -\T  - 1 50)3 
- 1. 26736  xl0~2(7"-150)2  +  12.2284(7" -150) +  5 1384] 
x  { (0.2409 tanh[9.5305  x  10  3(7")  - 12.259] }  +  0.7444  -D}21 
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3.3.3  Application  of  Sintering  kinetics  curves  on  the  fabrication  of  step-wise  FGMs  [191 


We  successfully  fabricated 

Zr02/50Zr02+50AI203/Al203  step-wise 
Functionally  Graded  Material  (FGM),  as 
shown  in  Fig. 20.  In  the  process,  the  powder 
selection  for  each  layer  is  based  on  the 
sintering  kinetics  comparison. 

To  achieve  the  flatness  in  the  FGM 
specimen  (no  significant  camber),  the 
sintering  kinetics  of  each  layer  should  be 
very  similar.  Comparing  Fig.  14  and  15,  one 
can  find  that  most  zirconia  powder  mixtures 
have  higher  densification  rates  than  the  Fi920-  A  representative  Backscattered 
alumina  powder  mixtures.  In  contrast  to  this  Electron  image  (BSEI)  of  the  three-layered 
trend,  the  curve  of  40%CERAC-2003+  Al203/Al203+Zr02/Zr02  FGM  [19], 
60%TZ3YS  (Powder  Z),  the  zirconia  powder  system  with  the  lowest  densification  rate, 
lies  between  100%CR-15  and  90%  CR-15+10%TMDAR  curves. 

Based  on  these  two  comparison  results, 
several  alumina  powder  mixtures  were 
tested  with  the  CR-15  powder  varying 
between  100%  and  85%.  The  TMA  testing 
results  revealed  that  the  alumina  powder 
mixtures  containing  94%  CR-15  and  6% 

TMDAR  (Powder  A)  has  very  similar 
densification  rate  curve  as  Powder  Z  (Fig. 

21 ).  Also,  sintered  samples  by  Powder  A  and 
Powder  Z  had  almost  the  same  final 
diameter.  These  two  powders  were  mixed  in 
equal  volume.  The  resulting  powder  mixture 
yielded  a  similar,  not  exactly  the  same, 
densification  curve  as  Powders  A  and  Z. 

This  should  be  due  to  the  change  in  the 


Fig. 21 .  The 
density  curves 
mixtures  (ramp 

particle  size  distribution.  However,  considering  the  excellent 
making  a  small  adjustment  in  the  TMDAR  proportion,  the 
CR-15,  4%  TMDAR,  and  20%  CERAC-2003,  30%  TZ3YS  by 
a  nearly  same  sintering  behavior  as  Powder  A  and  Powder  Z 
the  same  final  diameter. 


Relative  density 

densification  rate-relative 
of  the  three  chosen  powder 
rate:  10°C/min).  [19] 

sinterability  of  TMDAR,  by 
powder  system  with  46% 
volume  (Powder  M)  shows 
(Fig.  21)  as  well  as  almost 


The  shape  of  the  three  curves  in  Fig.  21  is  very  similar.  One  obvious  difference  is 
located  around  0.22-10'3/s  densification  rate  in  the  ascending  portion.  The  corresponding 
sintering  temperature  range  is  1150-1200°C.  There  is  a  plateau  for  layer  A  and  a  drop  in 
slope  for  layer  M,  but  there  is  no  noticeable  change  for  layer  Z  near  this  region.  As 
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explained  in  Section  4.3.1,  the  change  in  densification  rate  for  layers  A  and  M  is  a  result 
of  y-  to  a-phase  transition  of  alumina  [23]  since  1 0%  y-phase  alumina  is  present  in  CR-1 5. 
However,  since  layer  Z  does  not  contain  CR-1 5,  there  is  no  plateau  or  abrupt  slope 
change  in  its  densification  rate  curve.  The  similarity  among  the  sintering  kinetics  of  the 
three  powders  results  in  no  obvious  difference  in  shrinkage  rate  and  final  diameter 
among  layers  during  the  co-sintering  of  a  FGM  sample.  Therefore,  the  Powder  A,  Z.  and 
M  are  the  potential  powders  to  achieve  flat,  crack-free  three-layered  Al203/Zr02  FGM 
samples. 

3.3.4  CTEs  measurement  [29,  30] 

By  in  situ  reaction  of  WO3  to  Zr02,  we  fabricated  various  ZrW208/Zr02  composites,  as 
listed  in  Table  6. 

Table  6  The  W03/Zr02  mass  ratios  of  various  green  samples  and  the  corresponding 
resultant  ZrW208/Zr02  volume  ratios  in  the  sintered  samples. 


# 

W03/Zr02  mass  ratio  in 
reactant  powder 

ZrW208/Zr02  volume  ratios 
in  the  sintered  sample 

Final  relative 
density 

1 

0.159:1 

20:80 

77% 

2 

0.264:1 

30:70 

79% 

3 

0.38:1 

38.9:61.1 

80% 

4 

0.593:1 

51.5:48.5 

82% 

5 

1.096:1 

70:30 

83% 

6 

2.307:1 

90:10 

84% 

The  temperature-dependent  CTE 
value  of  each  sintered  sample  was 
measured  by  the  TMA,  and  shown  in  Fig. 
22.  As  shown,  both  Samples  1  and  2 
feature  positive  CTEs,  Samples  4-6 
exhibit  negative  CTEs,  and  the  CTE  of 
Sample  3  is  nearly  zero.  Fig.  22  shows 
that  the  CTE-temperature  curves  of  the 
composites  feature  a  shape  that  is  very 
similar  to  that  of  pure  ZrW208.  However, 
the  peak  caused  by  phase  transition 
occurs  at  a  lower  temperature  for  the 
composites.  This  “peak  shift”  is  possibly 
a  result  of  phase  transition  from 
a-ZrW208  (low  pressure  phase)  to 
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Fig. 22  Temperature  dependent  CTEs  of  the 
various  ZrW208/Zr02  composites  [29], 

y-ZrW208  (high  pressure  phase)  [31]  due  to  the  thermal  stress  induced  by  the  CTE 
mismatch  between  Zr02  and  ZrW208. 
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3.3.5.  Young’s  Modulus  measure  [29,  301 


The  TMA  unit  offers  a  facility  to  measure  45  >  - — - - 

Young’s  modulus  of  materials  using  the  4°- 
three-point  bending  method.  After  sintering  *  3-5- 

and  polishing,  the  samples  used  in  Young’s  S  3-0-  . 

modulus  testing  are  rectangular  ^  2.5-  T~ 1 — • — • — — r 

parallelepipeds  (cuboids)  with  edge  lengths  |  2.0- 
12.6mm,  4.2mm  and  2.5mm.  The  «  1 .5- 

O)  _ 

temperature  dependent  Young’s  moduli  of  5  10- 

the  fabricated  composites  listed  in  Table  6  ^  o,5-  -►-ZrWA 

as  well  as  that  of  ZrW208  and  Zr02  are  o.o-l — . — , — . — , — . — , — ■  ,  — , — . — , — . — 

,  _.  00  _.  _  .  ,,  0  100  200  300  400  500  600  700 

shown  in  Fig.  23.  Since  the  Zr02  in  the  T  4.  „ 

a  Testing  temperature  (  C) 

composites  is  processed  by  the  ramp/soak  Fig_23  Young’s  modulus  -  temperature 
path  with  quenching  step,  the  CERAC-2003  curves  for  composites  as  well  as  for 
for  Young’s  modulus  testing  was  also  ZrW208  and  Zr02  [29], 
partially  sintered  by  the  same  ramp/soak 

path.  Microcracks  formed  in  the  Zr02  samples  due  to  the  quenching  process.  The  highest 
testing  temperature  was  set  as  650°C  in  order  to  avoid  the  decomposition  of  ZrW208. 
The  room  temperature  Young’s  modulus  of  pure  ZrW208  was  found  to  be  4.31  GPa,  which 
nearly  matches  the  value  of  4.22GPa  reported  by  Chen  et  al.  [32],  However,  these  values 
quite  different  from  the  single  crystal  value  of  88.3GPa  reported  by  Drymiotis  et  al.  [22], 
Since  the  Zr02  sample  was  only  partially  sintered  and  contains  microcracks,  its  Young’s 
modulus  is  much  lower  than  that  of  a  fully-sintered  sample  (~200GPa)  and  it  does  not 
feature  a  large  decrease  from  100  to  400°C,  which  is  typical  for  fully-sintered  Zr02  [34], 
the  Young’s  moduli  of  the  composites  and  the  ZrW208  feature  a  value  change  around 
160°C  due  to  the  a  to  p  phase  transition. 


0  100  200  300  400  500  600  700 

Testing  temperature  (°C) 

Fig. 23  Young’s  modulus  -  temperature 
curves  for  composites  as  well  as  for 


3.3.6.  Verification  of  continuous  gradient  for  ZrO?-ZrW?Q8  continuous  FGMs  [35,361 


A  Zr02+W03  powder  mixture  and  Zr02 
powder  were  stacked,  co-compacted  and 
co-sintered  in  the  processing  steps 
commonly  used  for  multi-layer  materials. 
The  observation  of  the  cross-sectional 
microstructures  indicated  the  sintered 
samples  are  continuous  ZrW208/Zr02 
Functionally  Graded  Materials  (FGMs).  At 
the  same  time,  the  measurement  of  the 
radial  thermal  expansion  also  provided 
evidence  for  the  continuous  gradient  [35, 
36],  The  radial  surface  of  each  sample  was 
polished  to  achieve  two  parallel  planes.  The 
radial  thermal  expansions  of  each  sample 


-4-3-2-10  1  2  3  4 

_  Position  (mm) 

Fig. 24  The  radial  thermal  expansion- 

position  relationship  for  sintering  FGMs 
[35], 
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were  then  measured  at  several  locations  for  each  sample.  It  can  be  seen  in  Fig.  24  that 
the  radial  thermal  expansions  feature  a  continuously  changing  profile,  providing 
additional  evidence  of  a  continuous  gradient  of  the  sintered  FGMs. 
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